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parallel  feeds  are  presented.  To  obtain  the  formulas,  multiple  reflections  arc 
neglected,  and  devices  of  the  same  type  are  assumed  to  have  identical  electrical 
performance. 

The  approximate  results  were  compared  to  the  results  obtained  using  a 
scattering  matri.x  formulation.  Both  methods  were  in  agreement  in  predicting  RCS 
lobe  positions,  levels,  and  behavior  with  scanning.  I'he  advantages  of  the 
approximate  method  are  its  computational  efilciency  and  its  flexibility  in  handling 
an  arbitrary  number  of  coupler  levels. 
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I.  INTRODUCTION 


The  design  of  radar-steaithy  platforms  has  become  an  important  engineering 
problem.  The  principles  of  radar  stealth  have  been  well  known  for  several  decades, 
but  only  recently  have  technological  advances  allowed  practical  implementation  of 
these  principles.  Some  examples  of  stealthy  platforms  are  the  SR- 71  spy  plane,  the 
F-1 17A  fighter,  and  the  B-2  bomber.  All  of  them  have  low  radar  cross  section 
(RCS).  low  infrared  (IR)  emissions  by  control  of  heat  sources,  and  low  microwave 
emissions  by  using  "quiet"  radar  and  communications. 

The  RCS  of  future  military  platforms  will  be  lowered  significantly  by  the  use 
of  shaping  and  materials  selection.  Consequently,  attention  is  now  being  focused 
on  the  onboard  sensors  in  order  to  ensure  that  their  signatures  do  not  become 
predominant.  Of  particular  interest  are  wideband  phased  arrays,  which  can  possibly 
operate  at  the  same  frequencies  as  those  of  an  illuminating  threat  radar.  Thus  a 
high  performance  array  must  not  only  meet  the  system  antenna  operating 
requirements  (gain,  sidelohc  level,  etc.)  but  also  the  RCS  requirements.  It  is 
essential  that  an>'  technique  applied  to  minimize  the  radar  cross  section  does  not 
seriously  degrade  the  primary  operation  of  the  antenna  system. 


Solid  state  microelectronic  devices  now  permit  integrated  transmit/recei\e  or 
rcceivc-only  elements  to  be  collocated  with  array  apertures.  Self-calibrating  and 


adaptive  "smart  skin”  conformal  arrays  are  now  practical.  Increased  efficiency, 
reliability,  and  reduced  cost  are  the  benefits  of  this  technology.  This  technology 
makes  phased  arrays  more  appealing  than  other  antenna  types. 

In  characterizing  and.  subsequently,  minimizing  the  RCS  of  phased  arrays, 
reliable  computation  of  scattering  is  crucial.  It  is  necessary  to  consider  both  inband 
and  out-of-band  threat  frequencies  because  the  scattering  characteristics  of  an  array 
are  distinctly  different  in  the  two  frequency  regions.  If  the  incident  wave  is  in  the 
array's  operating  band,  the  radiating  elements  are  well  matched  and  threat  signal 
can  penetrate  into  the  feed  and  be  reflected  at  internal  mismatches  and  junctions. 
Even  for  a  well- matched  array,  there  can  he  a  large  number  of  scattering  sources 
that  add  constructively  under  some  conditions.  This  effect  depends  on  the  type  of 
feed  and  the  devices  incorporated  therein.  In  some  cases,  these  reflections  can 
significantly  modily  the  RCS. 

In  this  thesis,  the  inband  RCS  of  a  phased  array  w  ith  a  parallel  (corporate) 
feed,  as  shown  in  Figure  1.  is  examined.  The  main  objectives  arc: 

1.  to  find  approximate  equations  for  the  inband  RCS  of  phased  arrays  with 
parallel  feed  networks. 

2.  to  compare  the  approximate  solution  with  a  rigorous  solution  based  on 
scattering  parameters,  and 

3.  to  determine  the  RCS  behavior  for  various  feed  parameters  for  both  linear 
and  two-dimensional  arrays. 


In  general,  antenna  scattering  analysis  is  very  difficult.  Hansen  [Ref.  Ij 
describes  two  antenna  scattering  modes.  The  first  one  is  the  antenna  or  radiation 
mode,  which  is  determined  by  the  radiation  properties  of  the  antenna,  and  vanishes 
when  the  antenna  is  conjugate  matched  to  its  radiation  impedance.  1  he  second  is 
the  structural  mode,  which  is  generated  from  currents  induced  on  the  antenna 
surfaces.  The  two  modes  are  not  easily  identifiable,  in  particular,  when  the  array 
is  installed  on  a  platform.  In  this  thesis  only  the  antenna  mode  is  e.\amined.  which 
is  the  dominant  RCS  component  for  a  phased  array  in  its  operating  band. 

Chapter  11  provides  the  theoretical  background  for  array  antenna  RCS  analysis 
and  includes  the  derivation  of  the  approximate  formulas.  It  presents  information 
about  the  correspondence  between  RCS  lobes  and  the  location  of  periodic  scattering 
sources  within  the  feed  network.  Chapter  111  describes  the  formation  of  the  so- 
called  rigorous  solution.  This  method  includes  multiple  reflections  between  devices 
in  the  feed  w  hereas  the  appro.ximate  solution  only  considers  the  first  reflection  and 
neglects  higher  order  reflections.  Chapter  IV  deals  with  the  comparison  and 
analysis  of  results.  Finally.  Chapter  V  concludes  with  a  discussion  on  the  benefits, 
concerns,  and  recommendations  that  can  be  applied  to  the  analysis  of  low 
probability  of  intercept  radar  and  communications  systems. 


II.  THEORETICAL  BACKGROUND 


Antennas  with  identical  amplitude  and  phase  patterns  can  differ  in  the  way 
they  scatter.  The  possible  importance  of  this  point  in  connection  with  evaluation 
of  antennas  designed  to  have  identical  patterns  was  first  examined  by  R.  H.  Dicke 
I  Ref.  2],  who  proposed  "to  sec  what  can  be  done  in  the  way  of  differentiating 
between  a  good  and  bad  antenna  on  the  basis  of  scattering." 

A.  DEFINITION  OF  RCS 

The  assigning  of  a  radar  cross  section  to  a  gi\  cn  target  whether  it  is  an 
aircraft  or  a  ship,  is  based  on  the  fact  that  this  object  functions  as  an  antenna.  The 
back-scattering  cross  section  may  therefore  be  interpreted  as  a  measure  of  the 
antenna  currents  excited  on  such  an  object.  Radar  cross  section  is  a  measure  of 
power  scattered  in  a  given  direction  when  a  target  is  illuminated  by  an  incident 
wave.  (An  older  term  for  RCS  is  echo  area.)  Mathemaiicaib .  the  RCS  is  defined 


0(6.)  =  lim  4;tR^  — - - 

!£.(0,)P 


(1) 


where  R  is  the  distance  from  the  target  to  the  obser\'aiion  point  (receiver),  £,{9)  is 
the  scattered  electric  field  in  the  direction  of  the  receiv  er,  and  is  the  incident 
electric  field  (as.sumcd  to  be  a  plane  wave).  The  term  monostatic  means  that  the 
transmitter  and  receiver  are  co-localed  with  respect  to  the  target.  0  =  0,.  RCS  has 
units  of  square  meters. 

In  general.  RCS  is  a  function  of  the  angular  orientation  and  shape  of  the 
scattering  target,  as  well  as  frequency  and  polarization  of  the  transmitter  and 
receiver.  Some  typical  value.s  of  RCS  are  listed  in  Tabic  1. 

Table  1,  RCS  of  Common  Targets 


o  in  m- 

a  in  dBsm 

Target 

0,001 

-30 

Insects 

0.01 

-20 

Birds 

100 

20 

Fighter  Aircraft 

1000 

30 

Bomber  Aircraft 

10000 

40 

Ships 

6 


The  scattering  characteristics  tor  all  targets  fall  into  three  natural  regions. 
The  llrst  one  is  the  low- frequency  or  Rayleigh  region,  dell  net!  b>  .el.«  1  (1  is  the 
characteristic  length  of  the  target  and  k  -  2;t.'A  is  the  wave  niimberh  I  he  .second 
region  is  the  Mie  or  resonance  region,  where  il.  is  on  the  order  of  unity.  The  third 
is  the  high-frequency  or  optical  region,  where  ./\n  array  of  length  L  can 

fall  into  either  of  the  three  categories  depending  on  the  threat  radars  wavelength 
relativ  e  to  the  antenna’s  operating  wavelength. 

B.  SCATTERING  FUNDAMENTALS 

Antenna  scattering  has  been  the  subject  of  study  since  1950.  how  ever  there 
is  very  little  published  on  the  subject  in  the  open  literature.  Until  the  mid  I980‘s. 
work  has  concentrated  on  the  analysis  of  low  gain  antennas.  Recently,  high  gain 
antennas  have  received  attention  because  it  is  anticipated  that  they  will  be  in 
widespread  use  on  future  low  RCS  platforms. 

Tor  an  antenna  subjected  to  an  incident  electromagnetic  wave,  the  scattered 
field  has  traditionally  been  divided  into  two  modes;  the  antenna  mode  (£,)  and 
the  structural  mode  (f,)  [Ref  1  ].  The  total  scattered  field  (£)  is  the  sum 

£(Z^)  =  £  (ZJ  ^  £^  =  £/Z;)  >  It 


(2) 


vhere 


Z,  =  R^,  '  /.V„  is  the  radiation  impedance 
ll  =  effective  height  of  the  antenna 
Z,  =  3770 
K,  -  incident  field 

The  structural  mode  arises  from  currents  which  arc  induced  on  the  antenna 
and  the  surrounding  structure  when  the  terminating  load  is  equal  to  the  complex 
conjugate  t)f  the  antenna  impedance.  The  antenna  mode  results  when  the  induced 
current  delivered  to  the  antenna  feed  point  is  reflected  and  then  reradiated.  The 
antenna  mode  is  proportional  to  the  gain  of  the  antenna  for  a  given  direction  and 
a  modified  reflection  coefficient,  which  is  equal  to 

r  =  (3) 

z^-z; 

where  Z,  is  the  impedance  across  the  antenna  terminals. 

When  the  antenna  is  conjugate  matched,  the  modified  reflection  coefficient 
()-)  goes  to  zero  and  the  antenna  mode  vanishes.  Usual! \'.  this  is  achie\  ed  when  the 
antenna  impedance  has  a  real  value  and  the  terminals  are  connected  to  a  matched 
transmitter  and/or  receiver.  Consequently,  there  is  no  reflection  at  the  junction 


between  the  antenna  and  the  transmission  line.  However.  Z,  tor 


element  in  a 


array  is  dependent  on  the  angle  of  arrival  of  the  incident  wave  because  of  mutual 
ctiupling.  It  is  practically  impossible  to  force  r  to  zero  at  e\ery  angle 
simultaneously. 

If  there  is  a  mismatch  as  defined  by  the  modified  reflection  coelTicienl,  then 
antenna  mtide  reflections  occur.  The  amount  of  antenna  mode  scattering  depends 
on  the  gain  of  the  antenna  in  a  particular  direction.  Thus,  if  the  antenna  doe.s  not 
have  high  gain  outside  of  its  operating  band,  then  it  will  not  have  significant 
antenna  mode  scattering  even  though  tliere  is  a  large  modified  reflection 
coefficient. 

Kahn  and  Kurss  [Ref.  3]  iia\  e  demonstrated  that  for  a  large  class  of  antennas 
terminated  by  matched  receivers,  the  scattered  poyver  is  generally  greater  than  the 
absorbed  power,  equality  being  attained  for  minimum-scattering  antennas.  This 
result  has  frequently  been  interpreted  to  mean  that  no  conjugate- matched  antenna 
can  absorb  more  than  it  scatters.  This  implies  that  gain  must  be  lowered  to  reduce 
RCS.  However.  Green  [Ref.  4l  has  shown  that  an  antenna  can  ab.sorb  more  than 
it  scatters  if  the  gain  in  the  back  direction  exceeds  its  gain  in  the  forward  direction. 
This  is  the  case  for  all  high  performance  phased  atjay  antennas. 

C.  SCATTERING  CHAR.\CTER1STICS  OF  PHASED  ARRAY 

single  dipole  provides  low  directivity.  To  increase  antenna  .size,  and  hence 
directivity,  a  collection  of  elements  can  be  arranged  and  interconnected  to  form  an 
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arnsv.  tiic  Bisit  tLnneni  can  be  an  «i(>tntnt  (tf  iJtrf,  Immu,  fnicTxutriri  patchy  ipiftiJ, 
IT  (tipfi'le  rtcpasdin^  ml  iJle  Afijllriclllioti  la  ttn\a  prOvidl;  Vtn  illrttb^e 
Ilte  fjctdi  ttCtra  Ihr  array  cttmcirtS  must  uttaafai:  CDnslruciively  i,audj  Iti  Use  cttnusl 
dsremjons  mid  imEitbfj:  dcsvucEively  lancfik  eatli  fiiher>  in  Ilia  famaiHHie  ^fwice. 

rh*f  lUlli  of  EUi  iiiray  wuennlt  cm)  be  {bKJonJFWH-'d  into  ttic 
dCsjvtihs'U  in  lilt  [TCViMta  ^KJidii'  tbe  onEEnira  nLedc  anJ  Use  Stnietutral  ntudc.  fhc 
relHtivc  impartaniic  iif  Ihe  twj  tinns  will  lle|iC[|ii  pi  liiwsly  on  the  Ibrcm  ijcLjimiEy. 
As  siirtwi  in  Figure  -  Lhe  freiiuency  donsnin  is  separated  ipio  five  baiJs  For  a 
iveli’dcsipicd  aniflina.  lUc  RCS  h  Ihe  bund  sbfnild  be  low  becupse  mast 

uf  the  laTiJeal  caerE^  ti  deirvered  eo  Ibc  fiaUama  triad.  ]  LaweverT  even  LhiH^kS^h  Ihf 
imitvicUwsl  rejieefitiw  liotti  Tlic  juiechiui  flie  smallKii  array  caa  have  have  Sens 
(snlitHKiiflidf  oi  suelt  souteK.  rhm  ibt  KCii  eaa  itehrevc  signineaitT  levels  under 
some  eii>ncNliun&. 
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Figure  2,  Antcittsii  JVerjuefiey  naqti 


la  Usis  liitatb.  Cirtt^  Ihrnu  signals  in  the  ikpcnitfng  bflntl  of  the  jritenltU  wiH  He 
Lurvik JertO.  In  this  enat.  Hue  wave  ptnjctraici  pdIks  die  feed  netWrif^i  dlKl  t#  reJketed 
J|  ititemiU  iUilClinns  and  deviees.  Tbc  nmal  RCS  is  kfcferauacil  by  the  vector  illUi 


of"  the  individual  scattered  fields  that  return  to  the  aperture  and  rerad i ate.  These 
scattering  sources  arc  depicted  in  Figure  1  for  a  parallel  feed  network.  Tltc> 
include: 

1 .  aperture. 

2.  phase  shifter  inputs.  /  „ 

3.  inputs  of  the  first  level  of  couplers.  i\. 

•1.  loads  at  the  sum  and  difference  arms  of  the  first  level  of  couplers,  /"v. .  ''i. 
5.  sum  and  difference  arm  loads  for  higher  levels  of  couplers,  ctc. 

Typical  sources  of  scattering  are  the  following: 

1.  Devices  that  arc  not  perfect!}  matched  due  to  physical  limitations 
(manufacturing  and  assembly  errors,  variation  in  the  electrical  properties  of 
material,  etc.) 

2.  Surface  roughness  (assembly  tolerances,  discontinuities,  distortions,  etc.) 
j.  Errors  built  into  the  antenna  (e.s..  quantization  errors  from  subarrayina) 

4.  Edge  effecis- 

The  inband  antenna  mode  RCS  is  obtained  from  equations  (1)  and  (2)  and 


.summing  over  all  array  elements 


vhere 


r,(0.(p)  -  total  reflected  signal  returned  to  the  aperture  for  element  n 
when  the  wa\  c  is  incident  from  the  (O.ct)  direction 
k  ^  k{ujc  +  vy  tvi)  -  wave  vector 
u  =  sinGcostp 
V  =  sinGsimp 
W  =  COS0 

d„  -  position  vector  to  element  n 

=  effective  area  of  an  element  =  h'ZjAR^ 

=  normalized  element  scattering  pattern 

To  arrive  at  (4)  identical  elements  have  been  assumed;  the  variation  in  mutual 
coupling  near  the  array  edges  has  been  neglected.  This  allows  the  RCS  to  be 
separated  into  an  array  factor  and  an  element  factor,  just  as  in  the  radiation  case. 

Evaluating  (4)  requires  the  total  reflected  field  at  each  element.  A  rigorous 
solution  must  employ  a  network  matrix  formulation  such  as  scattering  parameters. 
If  multiple  reflections  within  the  feed  can  be  neglected,  an  approximate  solution 
can  be  obtained  by  tracing  signals  through  the  feed  and  hack  to  the  aperture. 


D.  approximate  method 


111  this  section  approximate  RCS  formulas  are  deriied.  The  following 
assumptions  are  made; 


1.  a!1  the  devices  of  the  same  type  are  assumed  to  have  identical  electrical 
characteristics.  That  means  that  all  the  radiaiina  elements  have  the  same 
reflection  coefficient  and  the  same  transmission  coefficient  None  of 
the  elements  is  ideal,  because  each  reflection  coefficient  r,  is  not  equal  to 
zero.  By  the  .same  token,  all  phase  shifters  have  a  reflection  coefficient  r.,. 
etc. 

2.  All  couplers  are  represented  by  magic  tees,  which  implies  equal  power 
splitting.  (This  is  not  a  low  sidclobe  feed). 

3.  In  the  operating  frequency  band,  all  feed  devices  are  well  matched  and 
therefore  higher  order  reflections  are  neglected  {r  «  1). 

4.  Only  scattering  from  the  aperture,  phase  shifter  inputs,  coupler  inputs,  and 
the  sum  and  difference  arms  of  the  first  and  second  levels  of  couplers  are 
considered.  Couplers  in  higher  levels  of  the  network  are  assumed  to  be 
perfectly  matched. 

5.  f.osslcss  dev  ices  arc  assumed  for  simplicit) .  which  implies 

r|*.|r|-  =  i  (5) 


for  a  device  where  r  is  the  reflection  coefficient,  and  t  is  the  transmission 
coefficient. 

6,  Identical  aperture  elements  with  a  Lambertian  scattering  pattern  (cos-0). 

7.  Edge  effects  are  not  included. 


8.  Assuming  that  only  one  scattering  source  dominates  at  any  given  angle,  the 
coherent  sum  of  the  scattered  signals  is  represented  b>  a  noncoherent  sum 


|£;  +  =  |£iP  *  \E^  -  -  ...  *  l£„  ^  (6) 


where  £„  is  the  reflected  signal  of  the  n-  element.  Thus,  the  total  RCS  can 
be  expressed  as  o  +  Oj.^  +  +  ... 

9.  Random  errors  are  neglected  since  they  only  contribute  to  an  average  RCS 
level. 


Parallel  feeds  are  suited  to  rectangular  clement  arrangements,  and  therefore 
linear  and  rectangular  array  geometries  will  be  studied.  Array  quantities  are 
defined  in  Figure  3.  Note  that: 


1.  For  the  case  of  linear  arrays,  all  elements  are  aligned  along  the  x-axis  and 
equalh  spaced,  d.  The  z-axis  is  broadside  to  the  array. 

2.  For  the  case  of  two  dimensional  planar  arrays  all  elements  arc  in  the  xy- 
planc,  uniformly  spaced  with  dimensions  d,.  and  d,.  and  numbers  of  elements 
N,.  and  N,  . 

3.  There  is  onl>  a  9  polarized  incident  field.  (For  linearly  polarized  elements 
in  the  xy  plane,  this  gives  rise  to  the  cos’0  scattering  pattern.) 

4.  The  phase  shift  per  element  introduced  by  the  phase  shifter  in  Figure  1  is 
A.  Furthermore,  the  phase  shifters  are  reciprocal. 
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0=1  n=^j 

Figure  3.  Two-Dimensional  Array  Geometr>' 


The  refleciion  sources  for  an  eight  element  linear  phased  array  with  a  parallel 
feed  are  shown  in  Figure  1.  The  incident  plane  wave  at  an  angle  0  arrives  at  each 
radiating  element,  which  is  the  first  scattering  source  encountered  with  rcllection 
coefficient  r,.  The  transmitted  signal,  which  is  determined  by  the  transmission 
coefficient  t,,  proceeds  to  the  phase  shifter.  Again,  if  the  phase  shifter  is  not 
matched  to  the  transmission  line,  a  reflected  signal  returns  to  the  apenure, 
depending  on  the  reflection  coefficient  .A  portion  of  this  reflected  field  is 
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reflected  again  by  the  radiating  eiement  since  it  is  assumed  to  be  reciprocal.  This 
is  a  second-order  reflection,  and  will  be  neglected  in  the  calculation  of  RCS.  The 
portion  of  the  signal  not  reflected  at  the  phase  shifter  is  transmitted  to  the  lirst  level 
of  couplers,  and  so  on. 

Assuming  that  each  phase  shifter  has  an  insenion  phase,  the  signal  that 
propagates  through  it  will  encounter  a  phase  shift  depending  on  the  antenna  beatn 
scan  angle  9^.  For  a  linear  phase  progression,  the  transmission  coefficient  for  the 
phase  shifter  at  element  u  is 

t  =  t  (7) 

Pn  P 

where  x„  =  ^'dsinO,  for  the  linear  array. 

Figure  4  shows  the  first  magic  tee  in  the  array  with  the  fourth  port  loaded. 
(It  can  represent  any  three-port  power  divider  as  well.)  Usually  mismatches  exist 
in  both  side  arms,  and  in  the  sum  and  difference  arms.  The  mismatches  result  in 
reflections  back  to  the  aperture,  with  reflection  coefficients  of  r.  at  ports  2  and  3. 
/V  at  pon  1,  and  r^  at  port  4.  The  angle  A  is  the  signal  phase  at  the  coupler  port 
3  relative  to  the  signal  phase  at  the  coupler  port  2.  and  includes  all  of  the  insenion 
phases  of  the  devices  between  the  couplers  of  the  first  level  and  the  apenure.  It 
also  includes  any  space  path  delay  reiative  to  the  origin  for  the  incident  wave.  If 
the  phases  of  all  reflection  coefficients  tc.xccpi  the  phase  shifter)  are  zero,  then 
=  a  -  Vm  where  a  =  tr/sine. 
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Figure  4.  Refiections  at  the  Ports  of  a  Magic  Tee. 

The  magic  tee  .shown  combines  signals  from  the  first  two  radiating  elements. 
Similarly,  the  next  magic  tee  combines  signals  from  the  third  and  founh  elements 
(2A  at  port  3  and  3A  at  pon  2).  The  portion  of  the  signal  that  enters  the  sum 
arm  propagates  down  to  the  second  level  of  couplers  where  some  reflection  occurs, 
with  the  remaining  signal  transmitted  to  the  third  level  of  couplers,  and  so  on.  For 
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a  parallel  feed  network,  the  number  of  radiating  elements  is  always  of  the  form  2"'‘ 
where  m  is  the  number  of  levels  of  couplers. 

The  signals  rellectcd  from  the  magic  tee  can  be  determined  from  its  scattering 
matri.\.  First  assume  that  the  magic  lee  is  perfectly  matched  and  therefore  has  the 
fo! lowing  scattering  matrix: 


0  110 

^  ^  1  0  0  1 

2  10  0-1 
0  1  -1  0 


(8) 


Referring  to  Figure  A.  the  input  signals  at  ports  2  and  3  are  E,  and  £, 
respectively.  Both  signals  have  unit  amplitude  but  differ  in  phase  by  A  radians. 
The  subscripts  on  E  refer  to  the  element  numbers  to  which  tee's  side  arms  arc 
connected.  In  the  sum  arm  the  combined  signal  is 


E,  1  .  2  2 

=  <j2e  ‘  cos(^) . 

Similarly,  for  the  difference  arm  the  combined  signal  is 

£  =  £  £,  =  il  .  yl 

i,  ‘^22 

-  y\/2«-''=  sin(  ^). 


(9) 


(10) 


18 


The  next  magic  tee  in  the  array  combines  the  input  signals  £.  and  from 
elements  3  and  4  Again,  both  signals  have  unit  amplitude  but  dilfer  in  phase  b\ 


A  (because  £j  has  a  phase  of  2 A  and  £,.  has  a  phase  of  3 A) 

=  E,  ^  E.  = 

I,  3  4  2  2 

■5  a 

—  A 

=  V  2  e  “  cos{ — )  » 

und 

E^  -  £3  £^  -  24  siii{y )  . 

In  general,  the  signals  in  the  sum  and  difference  arms  arc 

and 

E^  =  jy/2e  '  j 

for  n  =  1.2.  ....  N.'2  . 

The  reilecied  signals  from  the  sum  and  difference  arms  of 
couplers  have  the  form 


(12) 


(Oa) 


(13h) 


c  (irst  le\el  of 
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£r„  = 


(14a) 


and 

There  lore  ihe  total  reflected  signals  reiiimina  at  the  side  arm  inputs  are 


As  expected,  these  equations  indicate  that  the  RCS  depends  on  the  relative  phases 
of  the  signals  entering  the  side  arms. 

5  he  total  scattered  signal  due  to  first  level  of  eouplcrs  is  obtained  by 
summing  all  re Ilec lions  returned  to  the  aperture 
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(19) 


E.^  =  E[  ^  E.  >  E'i  * 


4- 


I  fit*  phase  shil'ter  deld>  and  path  del  as  in  the  direction  ol  tlie  observer  must  be 
included.  Summiiie  the  terms  in  (19)  yields  [Ref.  5  | 


■  sm{A/A)  ! 
^siii(2A)  j 


(20) 


This  can  be  interpreted  as  the  sum  of  returns  from  the  sum  arm  plus  returns  from 
the  difference  arms.  Applying  the  definition  of  RCS  gives 


4nA^  2  J  A 

- cos  — 

“  1 2 , 


I  /V 

^  sin(2A) 


(21a) 


sm(jVA)  - 
^sm(2A)| 


The  portion  of  the  si  anal  not  reflected  by  the  first  level  of  couplers  is 
transmitted  down  to  the  second  level  where  the  entire  process  is  repeated.  The 
total  signal  returned  to  the  aperture  due  to  reflections  at  the  second  level  of 
couplers  is 
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£/  =  4\/2cos^[— ]-[rj.cos^A  - 

UjL^:  4  J  sm(4A), 


(22) 


Again,  this  is  a  sum  of  RCS  contributions  from  the  sum  and  difference  arms. 


sinfAfA) 

a.2  (yj 

—  sm(4A) 
4 

47Ij4’  ^2  A  ) 

sin(//A) 

jsm(4A)| 

i'he  coupler  RCS  contribution  contains  three  factors.  The  first  one  is 
the  RCS  of  a  reflector  of  area  A  reduced  by  the  reflection  coefficient 
r.  The  squared  terms  in  large  brackets  are  an  array  factor  for  the  couplers. 
Final h',  the  remaining  factors  are  equivalent  to  an  element  factor  for  the  coupler 
sum  or  difference  arms. 

For  linearly  polarized  radiating  elements  along  the  \-axis  and  a  9-poiarized 
incident  wave,  the  element  factor  is 

=  cose .  (24) 

This  factor  can  be  lumped  with  A  to  form  a  projected  area.  The  phy  sical  area  is 


related  to  the  number  of  elements  and  the  spacing.  Thus  for  a  linear  arras 


A  -  Ndlco$B  (25) 

w  here 

A'  is  the  number  of  elements 
d  is  the  spacing  between  elements 
I  is  the  effcctise  height  ol'  the  element. 

For  every  RCS  contribution,  the  first  and  the  last  factor  of  f4)  remain  the  same,  but 
the  second  factor  is  different.  These  are  designated  as  array  factors  esen  though 
they  include  the  subarray  factors  for  the  couplers.  The  arras  factor  for  the 
radiation  element  contribution  is: 

'  'lf/sin(«)J 

The  array  factor  for  the  phase  shifter  contribution  is: 

^  2|si^]  (27) 

The  array  factor  for  the  side  arms  of  the  first  level  of  couplers  become,': 

.  2  7  jsin(//A)j  f28) 

liVsin(A)j 

fhe  arras'  I  actor  tor  the  sum  arms  of  the  first  level  of  coup  lets  becomes: 


Sm{NA) 

-^sin(2A) 


(29) 


which  occurs  in  (21a).  The  array  factor  for  the  difference  arms  of  the  fir.si  level 
of  couplers  becomes; 


siii(AfA)  ^ 
|sm<2A)J 


(30) 


which  also  occurs  in  (21b). 

The  array  factor  for  the  sum  arms  of  the  second  level  of  couplers  becomes: 


(31) 


which  is  related  to  (23a).  The  array  factor  for  the  difference  arms  of  the  second 
level  of  couplers  becomes: 

=  C^p^iCos2jA|sm‘(A) 
which  is  also  related  to  (23b).  By  combining  (4)  and  (6),  the  total  monosiatic 
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iiiband  antenna  mode  RCS  is 


o  =  (xf;  -  af';  ~af':  -  afI^  +  /if‘^  -af^^  af'd  .  (33) 

Consider  a  of  two-dimensional  planar  array  in  the  xy-piane  with  a  rcctanaular 
grid  as  shown  in  l  igure  3.  The  array  has  N\  by  N,  elements  spaced  by  d,, 
Parallel  feeds  are  used  to  combine  signals  for  elements  along  the  x-axis.  Thus  the 
results  for  the  linear  arras  can  be  applied  directly  to  the  two-dimensional  array  b\ 
choosing  A  properly  and  multiplying  all  terms  by  a  y-direction  array  factor.  If 
and  rp,  are  the  seanned  antenna  beam  angles,  then  let 

■  ^ii_jSm0jCos<p^  {34a) 

and 

Xoj.  =  A'dj.  sinG^sintp^  (34  b) 

be  the  intcrclement  phases  to  scan  the  antenna  beam.  Now  A  in  the  linear  array 
formulas  is  replaced  by  ^  where 

C.  =  -  Xox  •  <^5) 

All  scattering  terms  for  the  linear  array  must  be  multiplied  by  an  arras  factor  for 
the  s  -dimension.  For  scattering  sources  ahead  of  the  phase  shifters 


(36) 


l'sin(N,  W,,v)l 

^  f..  =  - - — ^ —  » 

[  jV^,sin(^:d^,v)  J 

and  for  scattering  sources  behind  the  phase  shifters 

[  A'3,sin(p  j 

where 


(37) 


Cy  =  .  (38) 

Thus,  for  a  two-dimensional  array  equations  (26)  through  (32)  become 


|'sin((V^/td^u)']  f 

'sin(A'3,;tri^.v)) 

’  [  A^^sin(/cd^M)J  [ 

2 1  sm{N^kd^uy 

1  [sinfA^^fed^.v)'] 

^  [A^_jSm(/td^u)^ 

H(V3,sin(A:d^,v) ) 

2^2^  (5in(^,QVsin(Wyg^ 

'  siii(^,y  ' 

—  sin(2  CP 


(sm(A^/p) 

1  A^.sin(g  J  ’ 


(39) 


(40) 


(41) 


(42) 
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smGV^C,) 

^sm{2C,) 


A^,sm(Q  1 


A  =  t~  t~  r^  cos‘[  1  coS"<Q 


C,'l  sin(/^,C,) 


-^sm(4Q 


'  sm(N^.  C,)'l 
,^,sm(Cp  ,  ’ 


1  ->  P’Crl  .  ’  siiiw  Q 

AF^^  =  cos‘  ~  sm'(Cp  - ^ 


-^sin(4  CP 


(sHNX)] 


Finall>,  ihc  total  RCS  is 


-  (AF;  ^AFp-  AF"  *  AF^  -  AFI  ^AF^^  AFIj  (46) 


where 


A 


.V ,  d^.  cose  . 


(47) 


III.  RCS  ANALYSIS  FOR  THE  RIGOROUS  SOLUTION 


In  the  approximate  method,  the  higher  order  reflections  are  neglected  because 
they  vary  as  r.  r-’",  etc.,  where  r  «  1.  A  rigorous  solution  based  on  a  scattering 
matrix  formulation  contains  the  effects  of  multiple  reflections.  Furthermore,  it  is 
possible  to  combine  scattering  matrices  with  the  method  of  moments  to  include 
interactions  between  the  feed  and  aperture  as  well  as  reflections  inside  the  feed, 
This  method  solves  the  problem  rigorously  by  obtaining  an  antenna  impedance 
matrix  that  describes  the  electrical  characteristics  of  the  antenna  surfaces  and  feed 
network.  The  method  of  moments  (MM)  impedance  matrix  is  combined  with  the 
feed  scattering  matrix  and  continuity  equations  that  relate  the  MM  expansion 
coefficients  to  the  feed  signals  as  described  in  [Ref,  6]. 

In  this  thesis  a  variation  of  the  rigorous  method  in  [Ref.  6]  is  used,  where  the 
radiating  element  is  represented  by  a  simple  two-port  device  with  reflection 
coefficient  /•,.  This  eliminates  the  method  of  moments  portion  of  the  antenna 
matrix  in  the  rigorous  solution  of  [Ref.  6],  and  the  problem  reduces  to  a  pure 
scattering  matrix  one.  This  allows  direct  comparison  of  the  rigorous  solution  with 
the  approximate  solution,  because  the  radiating  elements  are  modelled  the  same. 


A  typical  scattering  matrix  equivalent  network  with  eight  elements  is  shown 
in  I'igurc  5.  If  only  two  lc\  cls  of  couplers  are  modelled,  then  only  a  four-element 
array  need  be  considered.  Rejections  from  the  third  level  of  couplers  are  zero  and 
therefore  adjacent  four-element  subarrays  are  decoupled.  The  four-element  array 
problem  can  be  solved  by  the  scattering  matrix  formulation  in  [Ref  6|. 

In  the  rigorous  method,  the  number  of  radiating  elements  is  assumed  to  be  a 
multiple  of-l.  The  incident  plane  wave  field  incident  on  element  i  : 

5.  =  ,  (  =  1,  2,  „ ,  N  .  (48) 

The  N  elements  of  the  array  can  be  divided  into  N  '^  siibarrays  as  shown  in  Figure 
6,  where  the  a  coefficients  are  used  to  obtain  the  scattered  field  £' 

"  r  4  1 


liquation  (49)  has  been  obtained  using  the  fact  that  the  four-element  subarrass  are 
decoupled 


ft,  =  aj 

=  “6 
Cj  =  a- 

The  total  RCS  is  obtained  from  ( I ). 


(50) 
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4  -  element  linear  array 


A  A  A  *' 
f  f  f  V 

1  1  1 4 


Figure  6.  Scattering  Matrix  Equivalent  Network  when  the  Radiating 
Element  Match  is  Constant  with  Angle 


For  each  four-element  subarray,  there  are  eight  2-port  devices  (four  radiating 
elements  plus  four  phase  shifters)  that  yield  16  equations,  and  three  4-port  devices 
(t^vo  magic  lees  in  the  first  level  plus  one  magic  tee  in  the  second  level)  that  yield 
12  equations.  There  is  a  total  of  28  equations  and  28  unknowns  (28  o's)  for  the 
four  4 -element  subanray. 

The  imponant  characteristics  of  the  rigorous  method  are  the  follow  ing: 


1.  Interactions  between  all  devices  in  the  feed  and  radiating  elements  are 
included; 

2.  The  transmission  line  phase  between  devices  is  included.  Thus  the  effect 
of  changing  line  lengths  can  be  modelled. 
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IV.  RCS  OF  LINEAR  AND  TWO-DIMENSIONAL  ARRAYS 


In  order  lo  simplih'  the  analysis  and  reduce  the  demand  for  computer  time, 
only  linear  arrays  have  been  examined  using  the  rigonms  method.  Thus  a 
comparison  of  the  two  methods  is  only  given  for  linear  arrays.  RCS  contour  plots 
for  two-dimensional  arrays  were  also  obtained  to  illustrate  RCS  behavior  with 
beam  scanning.  The  MAIL  A  B  programs  are  included  in  the  Appendices. 

A.  RCS  DATA  FOR  THE  APPROXIMATE  SOLUTION 
1.  Linear  array 

For  the  approximate  method,  RCS  pattern  data  for  linear  arrays  was 
computed.  The  linear  array  program  computes  the  RCS  per  squared  wavelength 
in  dB  for  any  monostatic  angle  0,  number  of  radiating  elements  N,  and  scanned 
angle  0,,  The  contributions  from  each  scattering  source  can  be  broken  out 
individually  if  desired. 

I'he  inband  RCS  of  linear  arrays  with  16,  64.  and  128  elements  for 
9,  =  0  degrees  (no  scanning)  is  shown  in  Figures  7  through  9.  All  reflection 
coefficients  for  feed  devices  and  the  radiating  elements  are  equal  to  0.2 
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In  Band  RCS  of  a  linear  array  with  parallel  feed  N=  1 6  and  lhetas=0 


Figure  7.  Inband  RCS  of  a  Linear  Array  with  a  Parallel  Feed  for  N=16, 
e,  =  0°  (Approximate  Method), 
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Figure  8.  Inband  RCS  of  a  Linear  Array  with  a  Parallel  Feed  for  N=64. 
0j  =  O'’  (Approximate  Method). 


sigma/lamda'‘2  (dB) 


For  all  linear  array  calculations  the  spacing  is  =  Ul  and  the  effective  height  (or 
length  )  of  the  elements  in  the  v  direction  is  /  =  }J1  . 

Referring  to  the  figures,  note  that  the  number  of  major  lobes  (spikes)  is 
the  same  in  ail  cases.  For  a  small  number  of  elements  (e.g.,  N  =  16)  the  lobes  arc 
not  as  well  defined  because  they  are  broader  and  lower.  The  RCS  contributions 
from  first  and  second  levels  of  couplers  are  shown  in  Figures  10  and  11. 
respectively.  The  lobe  spacing  in  the  RCS  pattern  for  the  first  level  of  couplers  is 
determined  by  the  physical  spacing  of  the  couplers  {2d).  This  can  be  generalized 
for  higher  levels  also.  For  instance,  the  effective  spacing  of  the  second  level  of 
couplers  is  4rf.  Therefore,  as  more  levels  of  couplers  are  added  to  the  feed,  more 
lobes  appear  between  already  existing  lobes. 

Figure  12  illustrates  the  effect  of  beam  scanning.  Assuming  that  the  phase 
shifters  are  reciprocal  devices,  tlie  lobes  associated  with  mismatches  behind  the 
phase  shifters  scan  with  the  antenna  beam  because  of  the  factor  Xo  ■  The  large 
lobe,  at  0  =  45  degrees  in  Figure  12  is  due  to  in-phase  addition  of  the  scattered 
signals  passing  through  the  phase  shifters,  as  e.vpected.  Note  that  the  specular  lobe 
at  0  =  0  degrees  does  not  scan.  The  high  lobes  near  ±  85”  arc  due  to  Bragg 
diffraction  (the  RCS  equivalent  of  grating  lobes). 
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igTna/lamda''2  (dB) 


Figure  11.  Contribution  from  Second  Level  of  Couplers. 


in  Band  RCS  of  a  linear  array  with  parallel  feed  N=  64  and  thetas=45 


Figure  12.  inband  RCS  of  a  Linear  Scanning  Array  with  a  Parallel  Feed  for 
N  -  64.  ©5  =  45°  {Approximate  Method). 
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2.  Two  dimcmional  orray 

A  second  NlAlLAlt  pro^nun.  which  is  shown  in  Appendix  Ek  was  used 
III  generate  contnur  plots  of  inband  RCS  lor  IwiMtimensional  arrays  far 
imipllcify  .  only  sipinre  mrays  are  examined  l  N,  •  V,  luid  d,  •  U.  -  KC)  for  a 
spocifted  >CiUi  angle  (0,.  Aipiin  it  U  .is&umed  (bat  ail  reflection  cocflicuiitf^ 
arc  0.2  RCS  coiitoitrx  arc  shosVn  in  ligMrex  13  ihiough  16.  the)  arc  ploncd  jii 
dircciioD  cosine  space  (a  oud  v)  and  the  coniours  ericlo^  spatial  regions  of  iCS 
above  a  spccined  level  this  level  u  chosen  to  he  lO  OB  when  >V,  •  .V,  •  16  and 
20  dB  when  N,  *  N,  »  W  elonents 

rbk*  only  dillerence  between  the  arrays  ot  f  igures  13  and  1 4  ts  Ihc 
number  uf  elements:  16  by  16  versus  64  by  64  Both  pints  ore  symmetric  about 
the  horittonial  axis  at  v  0.  High  RCS  Is  present  along  the  principal  planes  of  the 
<UTa>'  and  then  drops  alt  away  from  these  axes  because  of  the  separable  pruiuci  in 
the  array  ruciors  1 1  can  be  seen  that  the  lobes  are  narrower  and  higher  as  the 
number  of  elements  itwreasc.  The  specular  lohe  appears  at  ihc  center  of  each  plot 
Figures  15  and  16  show  the  RCS  of  scanned  airnys.  The  fiest  tme  has 
H,  -  45  degrees  and  gi,  -  II  degrees,  and  tlie  second  B,  •  45  degrees  and  ig,  *  45 
degrees  Both  plots  are  symmetric  around  the  axis  v  =  0  As  in  the  linear  atray 
dasc.  the  lobes  onalnating  liom  rctloctions  behind  the  phase  shifters  «an  ainne 
wiU)  Uie  amoinu  beam 
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In  Band  RCS  of  a  two-dimensional  array  -  Contour  plot 


Figure  13.  Inband  RCS  of  a  Two-Dimensional  Array  with  a  Parallel  Feed  for 
N,  =  N,  =  16.  6,  =  0^  (Pi  =  0*.  Contour  Plot  at  10  dB  Level 
(.Approximate  Method). 
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in  Band  PCS  of  a  t^O'dimenaKonai  array  *  Contour  plot 


m  Band  RCS  of  atwo-dimensionai  array  -  Contour  plot 


In  Band  RCS  of  atwO'dimensional  array  -  Contour  plot 


Figure  16.  Inband  RCS  of  a  Two-Dimensional  Array  w'ith  a  Parallel  Feed  for 
=  Ny  =  64.  9j  =  45°,  (p^  =  45°.  Contour  Plot  at  20  dB  Level 
(Approximate  Method). 


B.  RCS  DATA  FOR  THE  RIGOROUS  SOLUTION 


I'or  the  rigorous  method,  pattern  plots  of  the  RCS  of  linear  arrays  are 
obtained  by  using  the  FOR  ("RAN  program  developed  in  (Ref.  6].  This  program 
computes  the  RCS  per  squared  wavelength  in  dB  for  any  monostatic  angle  0. 
Array  parameters  include  the  number  of  radiating  e!cment,s  A'  (A’  has  to  be  a 
multiple  of  8).  scan  angle  0,,  .  and  electrical  path  length  between  devices  1./^,  (in 
radians).  The  output  is  written  to  MATLAB  (lies  that  can  be  used  to  plot  the  RCS 
patterns. 

Figure  17  shows  the  RCS  values  of  a  linear  array  for  N  =  64  with  no 
scanning  (0^  =  0)  and  path  lengths  of  tpo  =  0  radians.  As  in  the  approximate 
method,  the  specular  lobe  at  0  =  0,  coupler  lobes,  and  Bragg  lobes  at  0  =  ±85 
degrees  are  evident.  The  same  array  is  examined  in  Figures  18  and  19  with  = 
t:.'4  and  respectively.  A  comparison  of  Figures  17  through  19  shows  that  the 
lobes  have  the  same  position  but  differ  in  amplitude.  SpecificalK.  for  \|v>  =  the 
specular  and  coupler  lobes  have  bigger  amplitude  than  for  ii/„  =  0.  However,  for 
ij/i,  =  ;r/4,  the  Bragg  lobes  have  smaller  amplitude  than  for  4/,,  =  0.  The  variation 
is  due  to  the  beating  of  mismatches  in  the  feed  The  addition  or  cancellation 
depends  on  the  line  lengths  connecting  the  devices,  q/,,.  This  effect  becomes  more 
complicated  as  the  beam  is  scanned  (Figure  20).  Some  lobes  ma\  disappear  at 
some  .scan  angles  because  the  phase  shifters  introduce  a  shift  that  cause  complete 
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sigma/lanida*2  {dB) 


sigma/lamda*2  (dB) 


in  Band  RCS  of  a  linear  array  with  parallel  feed  N-G4  and  thetas=45 


Figure  20.  Inband  RCS  of  a  Linear  Array  with  a  Parallel  Feed  for  N  =  64. 
9,  =  45°,  ut;  =  Tt/4  (Rigorous  Method). 
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cancellation  of  two  scattering  contributions.  This  is  not  predicted  by  the 
approximate  method  because  it  sums  the  individual  comribiiiions  non  coherently, 

C.  COMPARISON  SUMMARY 

The  rigorous  and  approximate  results  have  been  presented  for  broadside  and 
scanned  linear  arrays.  !n  both  cases  the  specular  lobes  have  almost  the  same 
inaanitude  (within  about  IdB)  but  the  coupler  lobes  vary  about  3  dB.  This  is 
attributed  to  the  noncoherent  summation  in  the  approximation.  When  the  beam  is 
scanned,  both  methods  predict  the  proper  lobe  locations,  unless  complete 
cancellation  occurs. 

From  a  practical  point  of  view,  the  approximate  results  are  very  close  to  the 
rigorous.  One  major  difference  is  the  computation  times.  For  N-  64.  the  rigorous 
method  takes  about  eight  times  longer  than  the  approximate  method  and  increases 
dramatically  when  .V  increases  (for  A'  =  128  over  12  times).  .A  more  important 
difference  is  that  if  the  number  of  levels  of  couplers  i.s  increased,  the  scattering 
equations  must  be  completely  rewritten  and  programmed-  For  the  approximate 
method,  more  terms  only  need  to  be  added  to  equation  (33). 
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V.  CONCLUSIONS 


Ail  approximate  scattering  model  for  arrays  with  parallel  feed  networks  has 
been  presented.  Calculations  for  several  eases  \\  ere  compared  to  a  rigorous  method 
which  includes  all  the  interactions  between  the  feed  devices  and  aperture.  The 
appro.ximate  method  was  in  good  agreement  with  the  rigorous  method  in  predicting 
RCS  lobe  positions,  heights,  and  behavior  with  scanning. 

There  are  several  advantages  to  the  appro.ximate  approach.  First,  it  is 
computationally  cOlcient.  allowing  two-dimensional  contours  to  be  generated  in 
minutes.  Second,  it  can  be  easily  extended  to  an  arbitrary  number  of  elements  and 
coupler  levels.  The  disadvantage  is  that  a  noncoherent  addition  of  terms  does  not 
predict  total  cancellation  conditions.  However,  most  RCS  designers  are  primarily 
concerned  with  ilie  "worst  ca.se"  conditions  for  highest  RCS.  and  in  this  sense  the 
approximate  method  is  sufficient. 

Future  efforts  should  be  directed  at  increasing  the  number  of  coupler  levels, 
and  adding  a  eoupling  network  in  the  ,v  dimension  of  the  two-dimensional  arrav . 
.Mso.  methods  of  reducing  the  inband  RCS  should  be  investigated. 
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APPENDIX  A 


MATLAB  PROGRAM  FOR  LINEAR  ARRAYS 

’/.  Phased  arrays 

•/.  Written  by  V.FLOKAS  -  22  APR,  1994 
X  Inband  RCS  of  a  linear  array  with  parallel  feed 

clg 

thetas=input(’ Enter  scanned  angle  in  degrees  thetas=  ’) 

X  Scanning  at  Phis=0 
this=thetas*pi/180; 
c=3eS ; 

f=input C ’ Enter  operating  frequency  in  Hz,  f=  ') 

lainda=c/f ; 

k=2*pi/lajiida; 

-<1=0 .  S*lajnda ; 
chi=k+d*sin(this) ; 
theta=linspace(-89 ,89,660) ; 
thi=theta*pi/l0O ; 

. alpha=k*d+sin(thi) ; 
j  6ta=alpha-chi«ones (1 , length (alpha) ) ; 

1=0 .5*lajiida; 

Ae=d*l: 

N=input( ’Enter  the  number  of  elements  N=  ’) 

A=N'»Ae; 
r=0 . 2 ; 

t=3qrt(l-(r"2)) ; 

al=(sin(N*alpha) ) ./(N*sin(alpha)) ; 
hl=f ind(isnan(al)) ; 
alChl)=ones(si2e(hl) ) ; 
a2=(sin(N*jeta)) . /(N»sin(jeta) ) ; 
h2=find(isnan(a2)) ; 
a2(h2)=one3(size(h2)); 
a3=(sin(N’‘jeta))  . /(  (N/2)*sin(2’f  jeta) ) ; 
h3=f ind(isnan(a3)) ; 
a3(h3)=onea(size(h3)); 
a4=(Btn<N*3eta))  ./((N/4)*sin(4t<jeta))  ; 
h4=f ind  C i sn  an ( a4) ) ; 
a4(h4)=ones(size(h4)) ; 

sigma=(4*pi*(A/lajnda)'2)*(<cosCthi))  .-2)  ,*((r-2)*(al.*2)  +  ((r*2)*(t*4) 
(al.-2))+((t'a)*(r'2)*(a2,"2))+((t'S)*(r*2)*((cos(jeta/2)).*4). 


Ca3.-2))+(Ct-8)*Cr-2)*((sin(jeta/2)).'4).*Ca3.-2))+ 
C(t~8)*(r''2)»((cos(jeta/2))  .  "4)  .  ♦  ((cos ( jeta) )  .  "4)  . -Ca4  "2))  + 
C(t''8)*(r~2)*C(cos(]Gta/2))  .  '4)  . * ((sin( jeta) )  .  *4)  ,*(a4.~2)))  ; 
sig=10*logl0(sigma/(lamda''2))  ; 
f igure(l) 

plot (theta, sig) ,grid 

xlabel( 'Monostatic  Angle  (deg)’) 

ylabGl(' sigma/ lajiida'2  (d8)  ' ) 

title(['In  Band  R.CS  of  a  linear  array  with  parallel  feed  H=  ’  ,num25tr(N)  ,  ’ 
and  thetas=' ,num2str (thetas)  ,  ]) 
gtext( 'Approximate  Method*) 
axis ([-90  90  -50  40]) 

'sec=(4+pi*(A/lamda)  “2)  ♦(  (cos  (thi)  )  .  *2)  .  *(  ( (t  “8)  ♦(r“2)*  (  (cos(jeta/2)  )  .  “4)  .* 

( (cos  (jeta)  )  .  "4)  .  *(a4.  '‘2)  )  +  (  (t“  3)  ♦(r~2)  »( (cos(  j  eta/2))  .  *4)  .»( (sin  (jeta)  )  .  "4)  . 
♦(a4.-2))); 

s=i0*logi0(sec/(lamda“2))  ; 
f igure(2) 

plot(theta,s)  ,grid.xlabel( 'Monostatic  Angle  (deg)  ')  ,ylabel( 'sigina/lamda'2 
(dB)') 

title(['In  Band  RCS  of  a  linear  array  with  Parallel  feed  M=  ’  ,n-um2str(N)  ,  ’ 
and  thetas=' ,iiura2str(theta3) ,  ]) 
gtext(* Second  level  of  couplers') 
axi3([-90  90  -SO  40]) 

f ir=(4*pi*(A/lamda) ■2)*((cos(thi) ) . *2) .♦( (t"8) «(r'2)*( (cos (j eta/2) ) . *4) . * 
((cos(jeta)).-4).*(a4.*2))+ 

((f8)*(r-2)*((cos(jeta/2))  .  *4)  .•((sin(jeta)  ."4)  .♦(a4.'2)))  ; 
sl=10*logl0(fir/(lajiida‘2))  ; 
figure (3) 

plot  (theta,  si)  ,grid,xlabel(  'Monostatic  Angle  (deg)  ’}  ,ylabel('sigTi(a/lainda'2 
(dB)-) 

title([’In  Band  RCS  of  a  linear  array  with  parallel  feed  N=  ’  ,num2str(N)  .  ' 
and  thetas= ' ,num2str(thetas) ,  ]) 
gtext( 'First  level  of  couplers') 
axis([-90  90  -50  40]) 
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MATLAB  PROGRAM  FOR  PLANAR  ARRAYS 

*/,  Phased  arrays 

•/.  Written  by  V.FLOKAS  -  22  APR.  1994 
X  Part  (a)  -  Linear  plot  for  2-d  array 
clear 
clg 

thetas=input( 'Enter  sceoined  angle  in  degrees  thetas=  ') 
this=thetas»pi/180; 

Phis=input( 'Enter  scanned  angle  in  degrees  Phis=  ') 

phis=Phis*pi/ 180 ; 

c=3e8; 

f=input( 'Enter  operating  frequency  in  Rz-  ’) 

iainda=c/f ; 

k=2*pi/lamda; 

dx=0.5*l3Lnida; 

dy=dx ; 

chiox=k*di*siii (this)*cos(phis)  ; 
chioy=k*dy*sin(this)*sin(phis) ; 
theta=linspace{-89,89.460) ; 
thi=theta*pi/180; 

PHI=Phis;  7.  IN  HONOSTATIC  RCS  (PHI  CUT  =  PHI  SCANNED) 

phi=PKlKpi/180; 

u=sin(thi) .*(ones(l .length(thi))*cos(phi))  ; 
v=sin(thi) .*(ones(l ,length(thi) )*sin(phi)) ; 
alpha=k'fdx»u; 
beta=k*dy*v ; 

jetax=alpha-chiox*ones(l ,length(alpha)) ; 
]etay=beta-chioy«ones(l .length (beta)) ; 
l=O.S»laiiida; 

Wx= input ( 'Enter  the  number  of  elements  Nx=  ’) 

Ny=Nx; 

A=Nx»Ny*dx*dy ; 
r=0,2; 

t=sqrt(l-(r'-2)); 

al=(sin(Nx*alp'na))  . / (Nx*sin(alpha) )  ; 
hl=f lnd(isnan(ai) ) ; 
al(hl)=ones(size(hl))  ; 
a2=(sin(Mx*jetax) ) . /(Nx*sin(ietax) ) ; 
h2=f  ind(isnaii<a2)  )  ; 
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a2Ch2)=ones(sizeCh2)): 

a3=(san(Nx*jetax) j . /( (Nx/2) *sin(2*j exax) ) ; 
h3=f  ir.d  {  i  snan  (  a3 ) )  ; 
a3(h3)=ones(size(h3)); 
a4=(sin(Ny*beta)) ./(Ny*sia(beta)) ; 
h4=f  jn(i(isnaii(a4)); 
a4{h4)=ones(size(h4))  : 
aS=(sin(Ny»j  etay))  .  /  (Ny*sin(j  etay)  )  , 
h5=find(isnaii(a5))  ; 
aSCh5)=ones(sizc(h5)); 

a6=(sin(Nx*jetax)) ./C(Nx/4)*sin(4*Tetax)) ; 
h6=find{isnan(a6)) ; 
a6(h6)=oties(size(h6)) ; 

sigma=(4*pi*(A/lamda)“2)*((cos(thi)) . ‘2) - *{ (r'2) « ( (al . *a4) . '2)+((r*2)+ 
(t"4)*((al.*a4).-2))+((t-8)*(r-2)*((a2.»a5).-2))+C(t-3)*(r-2)+ 
({ct)s(jetax/2))  .  *4)  .  * 

((a3.+a5)  .  "2) )  +  ( (t"8) ♦  Cr''2)  +  ((sin( jetax/2) )  .  ’4)  .»((a3.  »a5)  .  *2))  + 
((t'-8)*(r*2)*((cos(jetax/2))  .  *4)  .  *(  (cos( j etax)  >  . "4)  . ♦( {a6 .  ^aS)  . ‘2) )  + 
CCt* 8)*Cr-2)*((cos(jetax/2))  . ''4)  . «( etax) )  .'4)  .*({a6-*a5)  .  '2)))  ; 
sig=10+logl0  (sigma/  Clajnda'2)) ; 

'/.  Contribution  from  first  level  of  couplers  (sum  and  difference  arms) 
siginal=C4*pi+(A/lajnda)  “2)  +  ((co5(thi)) .  '’2)  .♦(((t'8)  +  (r"2)  + 

((cosCjetax/2)) .  "4)  . 

♦(Ca3.»a5).-2))+(Ct‘8)*<r-2)*((sinCjetax/2)) . *4) . » ( (a3 . taS) . ‘2) ) ) ; 
sigl=10*logl0Csiginal/{lamda'‘2))  ; 
y.  Contribution  from  second  level  of  couplers 

sigiiia2=(4+pi+(A/lamda)  ■2)*(  Ccos(thi))  .  *2)  .  »( ( (t'S) *(r'2)* ( (co3(jetax/2) )  . 
"4)  .*((cos(jetax)) ."4) .*((a6.*aS) . ■2))  +  ((t‘8)*(r'2)*((cos(jetax/2)) . "4) . 
*((sin(jetai)) . "4) . * ( (a6 . *35) . "2))) ; 
sig2=10*logi0(sigma2/(laiiida'2))  ; 
figure (1) 

plot (theta , s ig) , grid 
axi3([-90  90  -30  80]) 
xlabel( ’Honostatic  Angle  (deg)') 
yiabel(  '  sigina/lamda'2  (dB)  ' ) 

title(’ln  Band  RCS  of  a  planar  2-0  array  with  parallel  feed’) 
gtext(  [’Approximate  Method  -  Nx=' ,nm!i2str(Nx)  , '=Ny=N  and  thetas=’ . 
nuni2str(  thetas) 

gtext( [’Phis=' ,num2str(Phis) ,  ] ) 

pause 

figure (2) 


plot(theta,sigl) ,6rid 
axis ([-90  90  -30  80]) 
xlabel{’ Monostatic  Angle  (deg)') 
ylabel('sigma/lamda''2  (dB)  ' ) 

titleCln  Band  RCS  of  a  planar  2-D  array  with  parallel  feed  - 
First  Level  of  Couplers’) 

gt ext ( ['Approximate  Method  -  Nx=' ,num2str(Nx) , ’=Ny=N  and  t betas =’ 

, num2str (thetas ) ,  ]) 

gtext(  [’Phis=’ ,nuin2str(Phis)  ,  ]) 

figure (3) 

plot (theta, sig2) .grid 
axis {[-90  90  -30  80]) 
xlabel('Moiiostatic  Angle  (deg)') 
ylabel('3igraa/lanida''2  (dB)') 

titleCln  Band  RCS  of  a  planar  2-D  array  with  parallel  feed  - 
Second  Level  of  Couplers') 

gtext([’ Approximate  Method  -  ,num2str(Nx) , ’ =Ky=N  and  thetas= 

’ ,num2str (thetas) ,  ]) 

gtext(  ['Phis=' ,nuiii2str(Phis)  ,  ]) 

pause 

y.  Part(b)  -  Contour  plot  for  the  same  array 
clear  thi  phi  u  v 
du=0 . 005 ; 
dv=0.005; 
nu=fix(2/du)+l ; 
nv=f ix(2/dv)+l; 
for  i=l;nu 

u(i)*-i+{i-l)*du; 
for  j=l:nv 

v(j)=-l+(j-l)*dv; 

cc=(u(i))-2+(v(j))*2; 
cth(i, j)=l-((u(i)'2)+(v(j)-2)) ; 
if  (cOl) 

cth(i , j)=0; 

end 

[x,y] =raeshgrid(u,v) ; 
alpha=k*dx»x; 
beta=k*dy*y ; 
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jetax=alpha-chiox»ones{length(alpha) .length (alpha)) ; 
jeTay=beta-chioy*ones(length(beta) ,iength(beta)) ; 
ai=(sin(Nx*alpha)) ./(Nx*sin( alpha)) ; 

zl=f ind(isnan(al) ) ; 
al(zl)=ones(si2e(zl)) ; 
a2=(sin(Nx*jetax) ) , /(Nx*sln(jetax))  ; 

22=f ind(isnan(a2) ) ; 
a2Cz2)=ones(size(z2)) ; 
a3=(sin(Mx*jetax)) . / ( (Nx/2)*sin(2»] etax))  ; 
z3=f ind(isnan(a3) ) ; 
a3(z3)=ones(si2e(23)) ; 
a4=(siii(Ny»heta))  ./(Ny*8iii(beta)) ; 
z4=find(isnan(a4)) ; 
a4(24)=ones(size(z4) ) ; 
a5=Csir(Ny*jetay)) ./(Ny*sin(jetay))  ; 
z5=f  iad(isnaji(a5)  ) ; 
a5Cz5)=ones(size(z5))  ; 
a6=  Caiii(Nx»jet ax)  )  .  /  (  (Nx/4)  »sin(4* j etax) )  ; 

z6=f ind(isnan(a6)) ; 
a6C26)=ones(size(z6) )  ; 
q=C4*pi*(A/la]nda)  "2)  ♦  (cth)  ; 

n=(r-2)*((al.*a4).'2)+(r-2)*(t-4)*((al,*a4) .'2)+Ct*8)*(r-2)* 

((a2.*a5)  .'■2)  +  (t-8)»(r-2)*((cos(jetax/2))  -4).* 

{(a3.*a5).-2)+(t“a)*(r-2)*((sin(jetax/2)).-4).+((a3.*a5)--2)+ 
(t-8)*(r-2)*((cos(jetax/2)).“4).*((cos(jetax)).-4).*((a6.»a5).-2)  + 
Ct-8)*(r-2)*((cos(jetax/2))  .  *4) . * ( (sin(j etax))  .  *4)  .»((a6.=fa5)  .*2)  ; 
sigma=q. ♦n; 
sigma=abs (sigma) ; 

sig=10*logl0( (sigma/ (lamda'2))  +  eps»ones (length (sigma) .length (sigma)))  ; 

levs  [SO  20]  ; 

figure (4) 

aLxis(' square ' ) 

contour(sig.lev,u.v) .grid 

axis (’square ’ ) 

xlabel(’u’)  .ylabeK  'v' ) 

title('In  Band  RCS  of  a  two-dimensional  array  -  Contour  plot') 
gtextCC’ Approximate  Method  Nx= '  ,nuio2str (Nx) ,  '  =Ny=N  and  thetas®' 
,nuin2str (thetas)  ,] 

gtext(['Phis='  ,num2str(Phis) ,  ]) 
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